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We propose a way of detecting CP violation in a single neutrino oscillation channel at very long 
baselines (on the order of several thousands of kilometers), given precise knowledge of the smallest 
mass-squared difference. It is shown that CP violation can be characterized by a shift in L/E 
of the peak oscillation in the v e —Vp appearance channel, both in vacuum and in matter. In fact, 
matter effects enhance the shift at a fixed energy. We consider the case in which sub-GeV neutrinos 
are measured with varying baseline and also the case of a fixed baseline. For the varied baseline, 
accurate knowledge of the absolute neutrino flux would not be necessary; however, neutrinos must 
be distinguishable from antineutrinos. For the fixed baseline, it is shown that CP violation can be 
distinguished if the mixing angle #13 were known. 
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INTRODUCTION 

The vast majority of neutrino oscillation experiments 
0, 0, 0, 0, 0, la 0] can be understood within a three- 
neutrino framework @, 0, [l(| ■ Though oscillation phe- 
nomenology is becoming ever more quantitative, there 
remain unanswered fundamental questions regarding the 
nature of neutrinos. Oscillation studies will not be able 
to address some issues; that is, they cannot establish an 
absolute mass scale of neutrinos or whether neutrinos are 
their own antiparticles. However, oscillation experiments 
will be able to shed some light on the level of CP viola- 
tion in the lepton sector, the relative ordering of the mass 
eigenstates, and the octant of the $23 mixing angle. 

In the standard theory of neutrino oscillations, flavor 
interactions are not diagonal with respect to the mass 
eigenstates. Instead, flavor states are related to mass 
eigenstates via a unitary mixing matrix U which is typi- 
cally parametrized as in Ref. This parametrization 
contains three angles 6jk, which indicate the degree of 
mixing amongst the eigenstates, as well as three phases. 
The two Majorana phases cannot be probed in oscillation 
studies so that the Dirac phase S is of sole consequence in 
this realm. Finally, neutrino oscillations are not sensitive 
to the absolute mass scale but rather the mass-squared 
differences, Ajk = m 2 j—m\. Given three neutrinos, there 
are two independent mass-squared differences. A current 
analysis [HI indicates the one-cr values for the mixing 
angles 

6»i2 = 0.59±0.02, |6>i 3 |<0.09, 6» 23 = 0.76 ±0.08, (1) 
and the mass-squared differences 

A 2 i = 7.9±0.3xlCT 5 eV 2 , |A 31 | = 2.6±0.2 x l(T 3 eV 2 . 

(2) 

At the moment, there is no knowledge of the Dirac phase 
and only an upper bound exists on the magnitude of the 



mixing angle 6*13. Additionally, the lack of knowledge 
regarding the mass hierarchy is reflected in our ignorance 
of the algebraic sign of A31. 

Future experiments will more precisely determine these 
parameters as well as address our ignorance of the re- 
maining unanswered questions. Herein, we will focus 
upon the Dirac CP phase 8. CP violation can only be 
manifest in neutrino oscillations if there are at least three 
neutrinos with nondegenerate masses and nonzero mix- 
ing occurs among all three flavors. Given this, before 
searching for CP violation, one must first concretely es- 
tablish that 0i3 is nonzero. The CHOOZ reactor exper- 
iment [f| presently establishes the most stringent bound 
upon this mixing angle; future reactor experiments will 
achieve much greater sensitivities 12 , 13[ . Detecting CP 
violation will be difficult as its effects are modulated by 
the smallness of 613. As an added complication, for rela- 
tively long baselines through the earth, interactions with 
matter can result in differences between neutrino and an- 
tineutrino oscillations, even in the absence of intrinsic CP 
violation. This "fake" signal is not a consequence of a 
nontrivial CP phase but is due to the fact that the earth 
is made of matter (and not antimatter). Finally, igno- 
rance of the mass hierarchy can create difficulties when 
trying to extract intrinsic CP violation from such chan- 
nels. 

To ascertain the level of CP violation in the lepton 
sector, one might compare neutrino and antineutrino ap- 
pearance channels over the same baseline. This under- 
pins many experimental proposals; see Ref. [3] and the 
references contained therein. These experiments require 
long baselines; however, they are still in a regime in which 
the mass-squared dominance approximation is valid to 
first order in the ratio A21/A31. Should 6*13 not be too 
small, one could also, in principle, see the effects of CP 
violation within a single neutrino oscillation channel. In 
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fact, an index for characterizing the level of CP violation 
ascertained from a single-channel oscillation was devel- 
oped in Ref. [l^|, as an analog to the usual asymm etry 
index used when CP conjugate channels are operable [16j . 
Herein, we will consider the effect of CP violation upon a 
single oscillation channel. We begin by examining a neu- 
trino appearance channel in vacuum and then consider 
neutrinos which traverse a constant density mantle; the 
region of interest will require very long baselines in which 
terms involving A21 cannot be linearized. 



VACUUM OSCILLATION 

An understanding of the case of vacuum oscilla- 
tions will guide our efforts. We will use the standard 
parametrization of the PMNS mixing matrix U 



11]. The 



probability that an a-flavor neutrino of energy E will be 
detected as a /3-flavor neutrino at a baseline L is 

3 

V a(i {L/E) = Sap-i ]T 5R(C^)sin 2 (^ fe ) 

j>k 
3,fc=l 

+2 %( C ?k) M*Pik) » (3) 



j>k 
j , k = 1 



with C$ := U aj U* k Uf3 k U* 0J and <p Jk = l.27A jk L/E 
where the neutrino mass-squared differences Ajk = — 
m 2 . are in eV 2 and the ratio L/E is in units of km/GeV. 
Recalling that the CP phase changes sign when consider- 
ing antincutrinos, we see that the first sum in the oscil- 
lation probability consists of CP even terms whereas the 
last sum consists of CP odd terms provided a 7^ (3. In 
appearance channels a/ft the coefficients of each term 
in the sum are proportional to the Jarlskog invariant J 



m 



7/ 



With the standard parametrization one has 
J = S12C12S13C13S23C23S5 , 



(4) 



(5) 



where Sjk — sinOjk and so on. Given that the CP odd 
terms are modulated by S13, the mixing angle #13 must 
be appreciable to detect CP violation. If one were able 
to maximize terms linear in S13 then CP violating effects 
would be as large as possible. Note that to most eas- 
ily work with linear in #13 terms, we allow for negative 
#13 mixing angles as proposed in [ijj]- In Refs. IJj, 20( ] . 
we noted that such terms will be largest whenever the 
oscillation probability achieves a local maximum due to 
the smaller mass-squared difference A 2 i; as such, this is 
a fruitful region in which to explore the effects of CP 
violation. 



What then is the effect of the CP phase on a single 
appearance channel? For the CP even terms, the phase 
makes an appearance in the form of cos 5. Its size can 
modify the amplitude of the oscillation probability to or- 
der 6 2 . Of more interest are the CP odd terms. Sizable 
CP violation results in a phase shift of oscillation phases 
(pjk'i that is, for a given set of mixing angles and mass- 
squared differences, the peak oscillation probabilities will 
occur at a different value of L/E. If the phase shift is sig- 
nificant, then perhaps a broadband measurement around 
a local maximum would measure such a shift; or alter- 
natively, a varying baseline experiment could be used to 
determine the position of the maximum. This presumes 
sufficient previous knowledge of the mass-squared differ- 
ences, and this shift in the local maximum cannot inde- 
pendently resolve both #13 and sg- 

As the two mass-squared differences A21 ~ 8 x 1CP 5 
eV 2 and A 32 «2x 10~ 3 eV 2 are well separated, regions 
exist where oscillations can be approximated by a quasi- 
two-neutrino scenario. We first examine the region in 
which oscillations due to the A21 mass-squared difference 
is negligible. To first order, the remaining mass-squared 
differences are degenerate, A31 w A32. This near equal- 
ity causes the CP odd terms, given by, 



j>k 
3,^=1 



^(. C jk ) s[n ( 2L Pjk) ~ J[sin(2p 31 ) - sin(2<p 32 )] 



0. 



(6) 



to be zero to first order. 

A better region to explore is where the A21 driven 
oscillations are appreciable and the oscillations due to 
the larger mass-squared differences are unresolvable. In 
this region, we may approximate 

(sin 2 (^ 3j )> « i (sin(2^)) « 0, (7) 

for j = 1,2. Here, oscillations v e <-> are the dominant 
mode, so we will focus upon the appearance channel V e{i . 
To first order, this probability is 

7% w -25R(C 3 e f+C 3 e 2 1 )-45R(C7 2 e nsin 2 (^2i)-2Jsin(2^ 21 ). 

(8) 

The constant term due to the unresolvable higher fre- 
quency oscillations can be written as 



25ft(C 3 j' + C32 ) — 2s 13 c 13 s 2 3 



(9) 



this term is small as the mixing between mass eigenstates 
1 and 3 is small. The coefficient of the CP even term can 
be written as 

- 43?(C 2 e f) = 4[ S 2 2 c 2 2 c 2 3 ( C 2 3 - S 2 3S 2 3 ) + (c 2 2 - s\ 2 )K], 

(10) 

where we define K := si2Ci2Si3C 13 S23C23C«5- 



3 



We wish to characterize the effect of CP violation in 
a single oscillation channel V ejl . As such we will com- 
pare the consequences of maximal CP violation with no 
CP violation. First, let us determine the form of the 
oscillation probability whenever CP is conserved. We 
adopt the convention in Ref . ID, 21 1 in which we restrict 
the CP phase S € [0,7r) and allow negative mixing an- 
gle 6*i3 6 [— ir/2, 7r/2]. For no CP violation 5 = 0, one 
has J = and K = K max = si 2 ci 2 si3C? 3 s 23 c 2 3. Given 
this, with CP symmetry conserved, the Ve-v^ oscillation 
probability is 



A + B c pc sin 2 (^2i) 



(11) 



with A = 2s 2 13 c 2 13 sl 3 and B C pc = 4[s 2 2 c 2 2 c 2 3 (c 2 . 3 - 
s i3 s 23) + ( c i2 — s i2)^max]- Using the best fit mixing 
angles in Eq. ^ for #i 2 and O23 as well as the 1-a bound 
for 6»i3, one has A = 0.008 and B C pc = 0.47. As 9 13 is 
small, the contribution of the A" m ax term to Bqpc is on 
the order of 6%. The first oscillation maximum in this 
region occurs whenever ip2\ = ir/2. Using the best fit 
mass-squared difference in Eq. ([2|), one finds this maxi- 
mum at a baseline to energy ratio of L/E = 1.57 x 10 4 
km/GeV. The maximum value of the oscillation proba- 
bility is 0.48. 

When CP violation is maximal S = tt/2, one has J = 
Jmax = si2Ci2Si3cf 3 S23C23 and K — 0. Given this, with 



maximal CP violation, the oscillation probability is 
V e ^~ A + B sin 2 (ip 2 i ) + C sin(2</5 2 i ) 



(12) 



with A as above, B — 4s 2 2 c 2 2 c 2 3 (c 23 — s 2 3 s 23 ), and C = 
— 2J max . With some rearrangement, one may combine 
the two terms which are dependent upon the oscillation 
phase to acheive 



A' + B' sin 2 (v? 21 



(13) 



with the new terms 

A' = A + i(S - V-B 2 + 4C 2 ), B' = ^/B 2 +4C 2 , 

(14) 

and phase shift 



2C 

sin(20) = — . 



(15) 



Comparing with no CP violation, we note two dif- 
ferences: there is a slight shift in the maximum value 
of Vefj, and a shift in the value of L/E at which this 
maximum occurs. From the mixing angles in Eq. (fTj). 
one notes 4C 2 <C B 2 so that to a good approxima- 
tion A' « A — C 2 /B. Additionally, there is an adjust- 
ment in the depth of the oscillation probability; one has 
B' ps B + 2C 2 /B. Comparing the maximum value of the 
oscillation probability for the CP conserving (CPC) and 
maximally violating (CPV) cases, we find 



(p; 



CPC\ 



(Pg pv ) 



Acos(26 12 )K n 



C 2 
B ; 



(16) 



this difference represents a small shift in the maximum 
value of around 6%. 




L/E (10 km/GeV) 



FIG. 1: [color online] The curves represent v^-v^ oscillation 
probabilities in vacuum at a long baseline for the parameters 
indicated in the text. A detector resolution has been included 
in the curve. The (black) solid curve is the case in which CP 
is conserved. The (red) dotted curve is the case in which CP 
is maximally violated with #13 = 0.09. The (blue) dashed 
curve is the case in which CP is maximally violated with 
0i3 = -0.09. 

Perhaps more interesting is the value of L/E for which 
these maxima occur. This is found from the relative 
phase shift <f> between the two sine functions. Using the 
small angle approximation, one has 



C ( 2C 2 
B\ B 2 



(17) 



Keeping the first order terms in the mixing angle $13 and 
approximating maximal mixing for 623, one has 



713 



sin(20 12 )' 



(18) 



As expected, the shift in phase is proportional to the 
mixing angle #13. Using the usual mixing angles from 
Eq. ([1]), we find a phase shift of (j> = —5°. We note that 
should 0i3 = —0.09 (a negative angle), then the phase 
shift would be 4> — +5°. 

In Fig. [TJ we see that this phase shift results in a 
shift in the position of the oscillation peak. We com- 
pare the v e -v^ oscillation probabilities at very long base- 
lines for the CPC case and the two CPV cases. We 
use the mixing angles and mass-squared differences from 
Eqs. (P[5]) and mimic detector resolution so that the 
rapid oscillations due to the larger mass-squared differ- 
ence are present in the "wiggles" but are damped out. 
The peak probability for the CPC curve (solid line) oc- 
curs at L/E = 7r/(2. 54A 2i ) = 1.67 x 10 4 km/GeV as ex- 
pected. Relative to this, the CPV peak probability (dot- 
ted line) with positive 6*13 leads by an L/E of 0.18 x 10 4 
km/GeV, or lags (dashed line) by an equal amount for 



4 



negative 6*13. The height of the peak and the width of 
the peak decrese by a small amount for the CPV cases. 

In analyzing such data, it would be key to know the 
small mass-squared difference, A 2 i, with a high degree 
of precision. Additionally, we have taken for #13 its 1- 
a bound; if this mixing angle were smaller, the relative 
shifts in the two probability curves would be proportion- 
ally smaller and more difficult to detect. Clearly, to make 
measurements around L/E ~ 10 4 km/GeV one would 
need a very long baseline through the earth's mantle. For 
neutrinos of energy 0.5 GeV, for instance, this would oc- 
cur at a baseline of 5000 km. Matter effects then become 
quite significant so we next include them in the analysis. 



OSCILLATION IN MATTER 

Neutrinos which travel through sufficiently dense mat- 
ter undergo coherent forward scattering. Neutral current 
interactions produce no change in oscillation as all flavors 
interact equally. Charged current interactions, however, 
result in an effective potential for the electron flavor only 
[2H ]. In the flavor basis, the neutrino evolution equation 
in matter becomes 



id t Vf 



I 

2E 



UMU^ + V 



'7 



(19) 



where the mass-squared matrix is M, = diag(0, A 2 i, A 31 ) 
and the operator V operates on the electron flavor with 
a magnitude V = \/2GfN c , with Gf the Fermi coupling 
constant and N e the electron number density. For sim- 
plicity we shall assume matter of constant density. This 
will avoid the possibility of parametric resonances [23j 
which could obscure the features we wish to highlight. 
For a mantle density of 4.0 g/ cm 3 2J, |25| , the matter po- 
tential is about V ~ 1.5 X 10 -13 eV. For anti-neutrinos, 
one needs to change the algebraic sign of this potential. 
For anti-neutrinos, the MSW effect supresses the oscilla- 
tions and widens the peak so we will confine ourselves to 
the superior neutrino oscillation channel below. 

For producing results pictured in the graphs, we ex- 
actly solve Eq. (|T9|) . To best understand the origin of 
the effects seen, it is instructive to examine approximate 
analytical expressions for the oscillation probabilities. A 
constant density mantle is a relatively good approxima- 
tion; however, it limits our baseline to less than 5000 km 
[26| . This, in turn, places a limit upon the neutrino en- 
ergy if we are to reach the peak of the A21 oscillations. 
As a result, we will only consider neutrino energies below 
1 GeV. 

For such neutrinos traversing the earth, matter affects 
are most easily understood using the approximations de- 
veloped in Ref. 27 1; see also Ref. [2(J. With a constant 
density, one may diagonalize Eq. (fl~9| to achieve effective 
mixing angles in matter djf. as well as modified mass- 
squared differences A^,. For the current scenario, the 



mass-squared difference of most interest effectively be- 
comes 



A^ ps A 2 iy cos 2 26*12(1 — E/E R ) 2 + sin 26*12, (20) 



where the resonance energy is 

A 2 i cos 26*12 



Er = 



2Vcj 3 



(21) 



The remaining mass-squared differences are modified as 
well; however, the limit of incoherent oscillations used 
above for vacuum oscillations, Eq. ([7|), is still valid here. 
The mixing angle most severely affected by the matter is 
6*12; effectively, one has 



sin 20? 



sin 26*12 



cos 2 26 , 12 (1 - E/E R ) 2 + sin 2 26>i 2 



(22) 



The correction also results in a modified mixing angle 



'13 



'13 



sin26 ) i 3 £;i/ 

A^ 



(23) 



Finally, the remaining mixing angle #23 is unaffected by 
matter. We note that the only point at which a linear 
power of the mass-squared difference A31 enters is in the 
modification of 6*13. This correction for matter effects 
is quite small so that our results will be independent of 
the (unknown) neutrino hierarchy at the level of a 10% 
adjustment to the (already small) 6*13. We will assume 
normal hierarchy so that A31 > 0. 
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FIG. 2: [color online] Baseline for the first peak of the u^-v^ 
appearance probability driven by the effective mass-squared 
difference AJl. The resonant energy in the mantle is Er = 
0.10 GeV for the assumed neutrino parameters. The (black) 
solid curve is the case in which CP is conserved. The (red) 
dotted curve is the case in which CP is maximally violated 
with #13 = 0.09. The (blue) dashed curve is the case in which 
CP is maximally violated with #13 = —0.09. 

In this constant density mantle, the form of the os- 
cillation probabilities developed in Eqs. (jTT]-{T3J) is un- 
changed; one merely needs to replace all quantities with 
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their effective value in matter. From those expressions, 
we explore the consequences upon the very long baseline 
appearance channel. Given matter effects, the oscilla- 
tion phase ip 2 i now has additional energy dependence, 
thereby modifying the baselines needed to measure the 
peak of the V^ c curve. In Fig. we plot, as a func- 
tion of energy, the baseline needed to measure the first 
peak, ip Ca k = 7r£ | /(2.54 A21). In the plot, we express the 
energy in terms of the resonance energy; for the mixing 
angles and mass-squared differences assumed herein, the 
resonance energy in the mantle is around 0.10 GeV. For 
oscillations in vacuum, the curve would be linear; how- 
ever, in matter, the baseline L pea k decreases with energy, 
relative to vacuum oscillations, as the matter term begins 
to dominate the kinetic term in A^|. If CP were maxi- 
mally violated, then the baseline for the peak is shifted 
according to the phase shift <fi m which also carries an en- 
ergy dependence by virtue of the modified mixing angles. 
As with the vacuum case, a positive (negative) mixing an- 
gle 9i3 will result in a longer (shorter) baseline relative 
to the CPC conserved case; this is demonstrated by the 
dotted (dashed) curve in Fig. O 

The fact that the separation between CPC and CPV 
peak oscillation probabilities increases with energy is at- 
tributable to the energy dependence of the phase shift 
in matter <jf n . Improving on the 5° phase shift from 
the vacuum case would aid in resolving the size of the 
Dirac phase. We recall the definition of the phase shift, 
Eq. (fT5|) , and substitute the appropriate mixing angles 
in matter. The change in the mixing angles is most pro- 
found for 612 ; in fact, we see from Eq. (|22|) that this 
mixing angle vanishes as 1/E for large energies. This is 
seen in Fig. [3J 
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FIG. 3: [color online] The (black) solid line is the ratio 
sin sin 2&12. The (blue) dashed line is the phase shift 
(expressed in radians) in matter, <j) m , with #13 = —0.09. 

As a result, the term B m vanishes for large energies re- 
sulting in the limit <j) m — * n/A. This represents the max- 



imum phase shift between a CPC and maximally CPV 
peak. The dashed curve in Fig.[3]shows how 4> m increases 
for the energy range under consideration here. (Note that 
we use #13 = —0.09 so that we may plot a positive phase 
shift on the same graph.) The 612 factor contributes a 
term linear in energy which also combines with the cor- 
rection to #13 resulting in a small quadratic contribution. 
At E = 0.5 GeV, the phase shift is about 10°; for E = 1 
GeV, the shift has risen to 18°, thus enhancing the sen- 
sitivity to CP violation. 

In addition to the phase shift, the height of the peak 
value for the oscillation probability will be modified by 
matter effects. Recalling the coefficients Bcpc and B', 
we note that the amplitudes of the oscillations are dom- 
inated by sin 2 29^ in the energy region of interest. As 
such, the peak value for the appearance channel will de- 
crease with energy. With increased energies, matter ef- 
fects enhance the phase shift between the CPC and CPV 
cases; however, at the same time, the value of the oscil- 
lation maxima decreases, requiring a more intense beam 
or larger detector to maintain equal sensitivity. 



VARYING BASELINE 

Given the background of the approximate analytical 
framework above, we shall now show exact numerical 
results which demonstrate the features derived above. 
We continue with a constant density mantle and use the 
same oscillation parameters as above, but now we solve 
the propagation equation, Eq. (|19j) . sans approximation. 
Our interest is in the position of the peak oscillation in an 
appearance channel. Since locating a peak can be done 
without knowing the absolute normalization of the beam 
(typically a large systematic error in an experiment), this 
approach holds a distinct advantage. To locate a peak, 
a spread in L/E is required. This can be effected by 
varying the baseline for a fixed energy or measuring a 
broadband neutrino beam at a fixed baseline. 

We first consider a varying baseline for a narrowband 
neutrino beam. Sub-GeV atmospheric neutrinos would 
be the most likely source for such an experiment. Un- 
fortunately, a water Cherenkov counter such as Super- 
Kamiokande [H does not have the precision needed to de- 
termine the path length of lower energy neutrinos given 
the large scattering angle between the detected lepton 
and the incident neutrino. Additionally, such detectors 
do not have the ability to resolve neutrino from anti- 
neutrino which is here needed to achieve a clean signal. 
Future detectors, such as a magnetized iron calorimeter 
14j . might have sufficient ability to detect the sub-GeV 
atmospheric neutrinos while identifying whether they are 
particle or antiparticle. Despite the fact that the abso- 
lute flux of atmospheric neutrinos is only known to within 
15%, the proposed method of detecting CP violation is 
not sensitive to this as it is the shape of the oscillation 



L (10 km) 

FIG. 4: [color online] The curves represent v^-v^ oscillation 
probabilities through a constant density mantle at a varying 
baselines for neutrinos with a narrow energy band of E v = 0.2 
GeV. A detector resolution has been included in the curves. 
CP is conserved for the (black) solid curves. CP is maximally 
violated with #13 = 0.09 for the (red) dotted curves. CP is 
maximally violated with 6*13 = —0.09 for the (blue) dashed 
curves. 



probability curve which is the determining factor. Intro- 
ducing a finite detector resolution, we plot the appear- 
ance probability in Figs. 0HH] for three neutrino energies. 
We include the CPC case and the two CPV cases with 
#13 = ±0.09. Notice that the vertical and horizontal 
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FIG. 5: [color online] The same as Fig. [4] except E v = 0.4 
GeV. 

axes differ between the three curves. Upon examination, 
we note that the peaks of all curves are roughly located 
by the approximate baseline indicated in Fig. [21 the sepa- 
ration increases with energy. Note that the height of the 
peaks decreases with increasing energy, caused by the en- 
ergy dependence of sin 2 (26^); this decrease in amplitude 
effectively broadens the peaks. For energies above 0.5 
GeV, the peak of the CPC probability remains relatively 
near 4500 km while the peaks of the CPV curves move 



L (10 km) 

FIG. 6: [color online] The same as Fig. |4] except E v — 0.6 
GeV. 



farther away from this point. For instance, at 0.6 GeV, 
the peak of the CPV curve with 9 13 = -0.09 is 3900 km, 
and the peak for the CPV curve with 0x3 = 0.09 is around 
5100 km. We note too that there is a relatively signif- 
icant difference in the peak value of oscillation between 
the CPC and maximally CPV cases. How these features 
might be seen in an atmospheric neutrino experiment is 
left for future work. 



FIXED BASELINE 

Fixed very long baseline experiments have been previ- 
ously considered in the literature as a means to uncover 
CP violation, the value of 13 , and discerning the mass 
hierarchy [lj]. Typically, such experiments have base- 
line and neutrino energies that still permit a dominant 
mass-squared difference approximation to first or second 
order. As noted previously, we presently consider the 
other extreme in which the oscillations due to the larger 
mass-squared differences are unresolvable. As a source of 
neutrinos, we envision an accelerator beam stop. Using 
an additional detector near the source allows one to have 
a better estimate on the expected flux at a far detector. 
Also, one can potentially control whether the beam is 
running in a neutrino or antineutrino mode making it in- 
consequential as to whether the detector can distinguish 
the two. This fact would allow one to make use of exist- 
ing technologies, such as water Cherenkov detectors. In 
terms of our proposal, present technology would need to 
be extended to give lower energies in order to reach the 
needed values of L/E. 

To consider a fixed baseline detector, we will simulate 
a broad band neutrino beam with energies ranging from 
0.1 to 0.4 GeV with a flat spectrum. As we now are 
dealing with a range of energies, the effective mixing an- 
gles and mass-squared differences are no longer constant 
despite the constant density approximation. The phase 



7 



shift between the CPC and CPV cases will manifest it- 
self differently than in the fixed energy case. In Fig. [2l 
we plot the baseline for which the oscillation phase <y9 2 ™ 
achieves the value tt/2 in the constant density mantle, 
which simplified the analysis for fixed energy. There is 
no such simplification for a fixed baseline. 
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FIG. 7: [color online] The curves represent v^-v^ oscillation 
probabilities through a constant density mantle at a fixed 
baseline of 3500 km. CP is conserved for the (black) solid 
curve with Q\z = 0.09 and the (green) dot-dashed curve with 
#13 = —0.09. CP is maximally violated for the (red) dotted 
curve with d\s = 0.09 and the (blue) dashed curve with #13 = 
-0.09. 

The maximum value of the appearance channel not 
only depends upon the oscillatory phase but also 
upon the energy dependent amplitude of the oscillation. 
Recall, the energy dependence of the amplitude is domi- 
nated by sin 2 (26*^) which sharply decreases for energies 
greater than twice the resonant energy. On the other 
hand, for energies below twice the resonant energy, oscil- 
lation can be enhanced relative to the vacuum. This fac- 
tor must be included in determining the ideal baseline. 
Incorporating these two factors, we can determine the 
effect of baseline upon the detected neutrino spectrum. 
As an example, we examine a baseline around 4500 km. 
From Fig. [2J the oscillatory phases for the maximal CPV 
case with 6*13 = 0.09 and the CPC case achieve the value 
of tt/2 at 0.35 GeV and 0.62 GeV, respectively. However, 
as the CPC peak occurs at a higher energy its amplitude 
will be suppressed relative to the CPV case so that the 
two T ell curves could be readily distinguished. Turning 
to the CPV case with 6*13 = —0.09, we note that at this 
baseline the oscillatory phase ^21 is greater than 7r/2 for 
all energies under consideration. As a result, the oscilla- 
tion probability is further suppressed relative to the other 
two cases. In Figs. [?H9l we plot V e \i for three very long 
baselines. A horizontal shift (in L/E) of the location of 
the oscillation peaks is difficult to discern; however, there 
is a considerable relative vertical shift in the height of the 
curves. As this vertical shift is of import, we also plot 
the case in which CP is conserved but with negative mix- 
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FIG. 8: [color online] The same as Fig. [7] except with a base- 
fine of 4000 km. 
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FIG. 9: [color online] The same as Fig. [7] except with a base- 
fine of 4500 km. 



ing angle 6*13 = —0.09. In Ref. [19j, it was shown that a 
negative mixing angle will shift the oscillation probabil- 
ity vertically at these baselines and energies. Taken as 
a whole, there is no real signature amongst these curves 
which can charaterize both the level of CP violation and 
the sign of #13. However, if one assumed knowledge of 
6*i3, then there is a clear separation between the CPC 
and CPV curves. As an example, we consider the CPC 
and CPV curves for both having #13 = 0.09, the solid 
(black) curve and the dotted (red) curve at the baseline 
of 4000 km. The percent difference between the peaks of 
these curves is 14%. A similar relationship holds for the 
curves for 6*13 = —0.09, the dashed (blue) and the dash- 
dot (green) curves. Thus, should #13 become known and 
not be too small, this becomes a possible experiment for 
searching for CP violation. 



CONCLUSION 

We examined a means of detecting CP violation in a 
single neutrino oscillation channel. Such methods would 
require a non-zero value of #13 as well as a precise knowl- 



8 



edge of the small mass-squared difference. The CP phase 
was found to have the most profound impact for very long 
baselines at low energies, i.e., near the first peak of the 
oscillatory region of the smaller mass-squared difference. 
By examining a single oscillation channel, one avoids the 
possibility of a fake signal of CP violation attributable to 
matter effects. Additionally, for the energies and base- 
lines under consideration, mass hierarchy does not enter 
at a significant level. It was shown that CP violation 
can be characterized by a shift (in L/E) of the peak of 
the appearance channel relative to CP conservation, both 
in vacuum and with MSW matter effects. In fact, mat- 
ter effects enhance the shift at a fixed energy, and the 
effect is most pronounced for a varying baseline exper- 
iment rather than a fixed very long baseline. The best 
source for examining this effect is in atmospheric data, 
as fixing E and varying L is not so practical. For this 
accurate knowledge of the absolute neutrino flux is not 
necessary; however, neutrinos must be distinguishable 
from antineutrinos. Unfortunately, the matter-induced 
energy dependence of the parameters obscures the shift 
(in L/E) of the peak at fixed baselines. However, the 
height of such peaks exhibits significant dependence on 5 
so that the level of CP violation could be resolved if #13 
were known. Finally, the proposed measurements likely 
could not resolve the 613—6 degeneracy. 
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